Bose— Einstein condensation in trapped dipolar gases 
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We discuss Bose-Einstein condensation in a trapped gas of bosonic particles interacting dominantly 
via dipole-dipole forces. We find that in this case the mean-field interparticle interaction and, hence, 
the stability diagram are governed by the trapping geometry. Possible physical realisations include 
ultracold heteronuclear molecules, or atoms with laser induced electric dipole moments. 



Bose-Einstein condensation (BEC) of trapped atomic 
gases offers unique possibilities to highlight a gen- 
eral physical problem of how the nature and stability of 
a Bosc-condensed state is influenced by the character of 
interparticle interaction. In this respect, especially in- 
teresting are ultra-cold gases with attractive interaction 
between particles (scattering length a < 0). As known 
, spatially homogeneous condensates with a < are 
absolutely unstable with regard to local collapses. The 
presence of the trapping field changes the situation dras- 
tically. This has been revealed in the successful experi- 
ments at Rice || on BEC of magnetically trapped atomic 
7 Li (a = —14 A). As found in theoretical studies j|, 
if the number of Bose-condensed particles is sufficiently 
small (of order 10 3 in the conditions of the Rice experi- 
ments) and the spacing between the trap levels exceeds 
the mean-field interparticle interaction no\g\ ( n o is the 
condensate density, g = ATrTi 2 a/M, where M is the atom 
mass), there will be a metastable Bose-condensed state. 
In other words, the condensate is stabilized if the neg- 
ative pressure caused by the interparticle attraction is 
compensated by the quantum pressure imposed by the 
trapping potential. In some sense, this is similar to the 
gas-liquid phase transition in a classical system with in- 
terparticle attraction: The gas phase is stable as long as 
the thermal pressure exceeds the (negative) interaction- 
induced pressure (see j|). 

The recent success in creating ultra-cold molecu- 
lar clouds opens fascinating prospects to achieve 
quantum degeneracy in trapped gases of heteronuclear 
molecules. In a sufficiently high electric field "freezing" 
their rotational motion, these molecules interact via the 
dipole-dipole forces. This interaction is long-range and 
anisotropic (partially attractive), and there is a non- 
trivial question of achieving BEC and manipulating con- 
densates in trapped gases of dipolar particles. 

Thus far, only the interaction between (small) atomic 
dipoles has been included in the discussion of the con- 
densate properties. Goral et al. || considered the ef- 
fect of magnetic dipole interaction in a trapped spin- 
polarized atomic condensate. Magnetic dipoles are small 
(of the order the Bohr magneton /ib), and even for atoms 
like Chromium (6/ib) the magnetic interactions are dom- 



inated by the Van der Waals forces. Nevertheless, for a 
relatively small scattering length a the condensate wave 
function may develop novel structures reflecting the in- 
terplay between the two types of forces. These effects 
can be amplified by modyfing a, which hopefully will 
soon become a standard technique P], and could even- 
tually appear in other systems, such as polar molecules, 
as pointed out in Ref. B. Similar effects have been dis- 
cussed by Yi and You |n| for ground-state atoms with 
electric dipole moments induced by a high dc field (of the 
order of 10 6 V/cm). These authors have demonstrated 
the validity of the Gross-Pitaevskii equation (CPE) for 
this system, constructed the corresponding pseudopoten- 
tial, and determined an effective scattering length. 

In this Letter we discuss BEC in a trapped gas of dipo- 
lar particles, where the interparticle interaction is dom- 
inated by the dipole-dipole forces. Possible realizations 
include the (electrically polarized) gas of heteronuclear 
molecules as they have large permanent electric dipoles. 
We also propose a method of creating a polarized atomic 
dipolar gas by laser coupling of the atomic ground state 
to an electrically polarized Rydberg state. Similarly to 
the condensates with a < 0, dipolar condensates are un- 
stable in the spatially homogeneous case, and can be sta- 
bilized by confinement in a trap. However, we find a 
striking difference from common atomic condensates: In 
the BEC regime the sign and the value of the dipole- 
dipole interaction energy in the system is strongly in- 
fluenced by the trapping geometry and, hence, the sta- 
bility diagram depends crucially on the trap anisotropy. 
This offers new possibilities for controlling and engineer- 
ing macroscopic quantum states. Remarkably, for dipoles 
oriented along the axis of a cylindrical trap we have 
found a critical value l* — 0.4 for the ratio of the ra- 
dial to axial frequency I — (ujp/uJz) 1 ^ 2 : Pancake traps 
with I < I* mostly provide a repulsive mean field of the 
dipole-dipole interaction, and thus the dipolar conden- 
sate in these traps will be stable at any number of parti- 
cles N. For I > Z» the stability of the condensate requires 
N < N c , where the critical value N c at which the collapse 
occurs is determined by the condition that (on average) 
the mean-field interaction is attractive and close to ui p . 

We consider a condensate of dipolar particles in a cylin- 
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drical harmonic trap. All dipoles are assumed to be 
oriented along the trap axis. Accordingly, the dipole- 
dipole interaction potential between two dipoles is given 
by V d (R) = (d 2 /R 3 ){l - 3cos 2 6>), where d is the dipole 
moment, R the distance between the dipoles, and the 
angle between the vector R and the dipole axis. The 
dipole-dipole interaction is long-range, and one can no 
longer use the pseudopotential approximation for the 
mean field. Similar to |Q, we describe the dynam- 
ics of the condensate wave function ?/>(r, t) by using the 
time-dependent GPE 

+ s »(r,nr + d 2 J <n" 1 ~ 3 _ c ^ 3 e \t(i",t)\'\mty (i) 

Here ^(r, t) is normalized to the total number of conden- 
sate particles N. The third term on the rhs corresponds 
to the mean-field interaction due to (short-range) Van 
der Waals forces, and the last term to the mean field of 
the dipole-dipole interaction. Assuming that the inter- 
particle interaction is mostly related to the dipole-dipole 
forces and d 2 ^> \g\ = ATrh 2 \a\/M , we omit the (Van der 
Waals) term g\tp(r, t)\ 2 ijj(r, t). 

The wave function of the relative motion of a pair of 
dipoles is influenced by the dipole-dipole interaction at 
interparticle distances |r — r'| < = Md 2 /h 2 . This 
influence is ignored in the dipole-dipole term of Eq. (0) , 
as the main contribution to the integral comes from dis- 
tances |r — r'| of order the spatial size of the condensate, 
which we assume to be much larger than r* . 

As mentioned above, in the spatially homogeneous case 
the dipolar condensate is unstable. For all dipoles paral- 
lel to each other, by using the Bogolyubov method one 
easily finds the anisotropic dispersion law for elementary 
excitations: e(k) = [E% + 8nE k n d 2 (l/3 - cos 2 9k)} 1 / 2 , 
where E k = h 2 k 2 /2M, no is the condensate density, and 
Ok the angle between the excitation momentum k and 
the direction of the dipoles. The instability is clearly 
seen from the fact that at small k and cos 2 Ok > 1/3 one 
has imaginary excitation energies e. 

To understand the influence of the trapping field on 
the behavior of the dipolar condensate, we have numeri- 
cally simulated Eq. (Jl|) for various values of the number of 
particles N, dipole moment d, and the trap aspect ratio 
I. By evolving Eq.(|l|) in imaginary time, we have found 
the condition under which the condensate is stabilized 
by the trapping field and investigated static properties 
of this (metastable) Bose-condensed state. 

For the stationary condensate the wave function 
ip(r,t) = ipo(r) exp (— i fit /ft), where u is the chemi- 
cal potential, and the lhs of Eq.(Q) becomes fxipoi^)- 
The important energy scales of the problem are the 
trap frequencies oj z , u> p and the dipole-dipole interac- 
tion energy per particle, defined as V = (1/N) J Vd(r — 
r')^ 2 (r)-0 2 (r')drdr'. Accordingly, the quantity V/Tiu> p , 
the aspect ratio of the trap I, and the (renormalized) 



number of particles a = JVr,/a max , (with a max = 
(h/2Mu> m i n ) 1 / 2 being the maximal oscillator length of 
the trap) form the necessary set of parameters allowing 
us to determine the chemical potential and give a full de- 
scription of the behavior of a trapped dipolar condensate. 
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FIG. 1. Critical value a c = A^ c r t /o max and (in the inset) 
the corresponding condensate aspect ratio L c versus I. 

We have found that the dipolar condensate is stable 
either at V > 0, or at V < with \V\ < fiuo p . This re- 
quires N < N c , where the critical number N c depends on 
the trap aspect ratio I. The calculated dependence N c (l) 
is presented in Fig. [l] and clearly indicates the presence 
of a critical point Z* = 0.4. In pancake traps with I < 1* 
the condensate is stable at any N, because V always re- 
mains positive (see Fig. ||). For small N the shape of 
the cloud is Gaussian in all directions. With increasing 
N, the quantity V increases and the cloud first becomes 
Thomas-Fermi in the radial direction and then, for a very 
large N, also axially. The ratio of the axial to radial size 
of the cloud, L = L z / L p , continuously decreases with in- 
creasing number of particles and reaches a limiting value 
at iV — ► oo (see Fig. ||). In this respect, for a very large 
N we have a pancake Thomas-Fermi condensate. 

For I > 1 the mean-field dipole-dipole interaction is 
always attractive. The quantity \V\ increases with iV" 
and the shape of the cloud changes (see Fig. || and Fig. 
||) . In spherical traps the cloud becomes more elongated 
in the axial direction and near N = N c the shape of the 
cloud is close to Gaussian, with the aspect ratio L = 2.1. 
In cigar-shaped traps (! > 1) especially interesting is 
the regime where hu z <C \V\ <C hui p . In this case the 
radial shape of the cloud remains the same Gaussian as 
in a non-interacting gas, but the axial behavior of the 
condensate will be governed by the dipole-dipole inter- 
action which acquires a quasilD character. Thus, one 
has a (quasi) ID gas with attractive interparticle inter- 
action and is dealing with a stable (bright) soliton-like 
condensate where attractive forces are compensated by 
the kinetic energy (TTJ . With increasing TV, L z decreases. 
Near N = N c , where \V\ is close to huj p , the axial shape 
of the cloud also becomes Gaussian and the aspect ratio 
takes the value L w 3.0. For I* < I < 1, the dipole- 
dipole interaction energy is positive for small number of 
particles and increases with N. The quantity V reaches 



2 



its maximum, and the further increase in N reduces V 
and makes the cloud less pancake. At the critical point 
N = N c the shape of the cloud is close to Gaussian and 
the aspect ratio L < 3.0, tending to 1 as I — > Z*. 
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FIG. 2. Dipole-dipole interaction energy V versus a for (a) 
I = 10, (b) 2 = 1, and (c) I — 0.1. In the inset the figure (c) 
is depicted in a larger scale. 



where B(L) = 2 + L? - 3L arctan[Vl - i 2 /i]/\/l - L 2 . 
Similarly, one can find the corresponding expressions for 
a, L z and L p as a function of ui z<p . The result of Eq.(j^) 
differs by less than 15% from our numerical calculations. 
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FIG. 4. Condensate widths (in units of a z ) in the dipole 
direction (solid line) and perpendicular to the dipole (dashed 
line), for the collapse in a spherical trap of frequency w 
(ct(0) = 11.5 < cr c = 11.8, and a(t) = 11.5 + 0.8wt). 
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FIG. 3. Condensate aspect ratio L versus a for (a) I — 10, 
(b) I = 1, and (in the inset) (c) I — 0.1. Figures (a) and (b) 
end when the system enters the unstable regime. 

To gain insight in the nature of the stability of dipolar 
condensates in pancake traps we performed a variational 
ansatz assuming a Gaussian shape of the cloud (c.f. [[Tol): 



V>Q = W 1/2 (2^)- 3 / 4 (L^)-l/2 e -p 2 /4^ e -, 2 /4L?. 
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Minimizing the energy functional H of the system, we 
found the aspect ratio of the cloud L and established that 
it decreases with increasing a for I < Z», and increases 
otherwise. The point Z* can be estimated by requiring 
dL/da\ a= o — 0, which provides the value Z» = 0.41 in 
good agreement with the numerical calculation. 

For understanding the behavior of the dipolar conden- 
sate near the critical point N — N c , we observe that 
at this point the local minimum of H becomes a saddle 
point. Hence, at N c one has (d 2 H/dL 2 z )(d 2 H/dL 2 p ) - 
(d 2 H/dL z dL p ) 2 = 0, in addition to 8H/dL z = 
dH/dLp = 0. This gives the relation between L and 
I at the criticality: 



(2L 2 + 1)(5 + 10Z 4 ) 6B(L){1 + 21 4 



2L 4 + 1 



(L 2 - l) 2 



-1 = 2- 



L 2 ' 1 



(3) 



We have also analyzed the dynamics of the instability, 
evolving (in real time) the initially stable condensate by 
a slow increase of a. In our approach based on the GPE, 
similarly to the case of a < (see [^2[ ), the dipolar con- 
densate collapses to a point on a finite time scale. The 
ratio L increases moderately in the course of the collapse. 
The results for the spherical trap are presented in Fig. 
[|. The collapse in cigar-shaped traps occurs in a similar 
way, since the initial shape of the collapsing dipolar cloud 
is almost the same as in spherical traps (see Fig.l). 

As we see, the ground state of a dipolar gas exhibits a 
very rich behavior and one finds various BEC regimes. In 
order to electrically polarize an ultra-cold cloud of het- 
eronuclear molecules, and thus create a molecular dipo- 
lar gas one should have the electric field which provides 
the Stark potential dE greatly exceeding the spacing be- 
tween the lowest rotational levels of the molecule. Then 
the rotational motion of the molecules will be "frozen" 
and their dipole moments will be oriented along the di- 
rection of the field. For most of the diatomic molecules 
the rotational level spacing is in the range from 0.1 to 1 
K, and hence the required electric field is E ~ 10 3 V/cm. 

At present, molecular condensates are not achieved ex- 
perimentally. We thus propose and analyze an alterna- 
tive method of inducing electric dipole moments, which 
can be used in atomic condensates. The idea is to apply a 
constant electric field and to optically admix the perma- 
nent dipole moment of a low-lying Rydberg state to the 
atomic ground state. Rydberg states of hydrogen and al- 
kali atoms exhibit a linear Stark effect : in hydrogen, 
for example, an electric field E s splits the manifold of Ry- 
dberg states with given principal quantum number n and 
magnetic quantum number m into 2(n — \m\ — 1) Stark 
states. The outermost Stark states have (large) perma- 
nent dipole moments dp> ~ n 2 eaB (with clb the Bohr ra- 
dius), and there will be an associated dipole-dipole force 
between the atoms. 
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This dipole-dipole interaction can be controlled with 
a laser [Q . This is achieved either by admixing the per- 
manent dipole moment of the Stark states to the atomic 
ground state with an off-resonant cw laser, or by a stro- 
boscope excitation with a sequence of laser pulses. The 
pulses tuned to the lowest Stark state of a given Ryd- 
berg manifold should be separated by the time T, have 
duration 2 At < T and area 2ir @. The field E s and 
the "dressing" light have to be chosen such that they do 
not couple the selected lowest Stark state to other states, 
and the spacing (~ nea^Eg) between the adjacent Stark 
states should greatly exceed the mean-field dipole-dipole 
interaction in order to avoid the interaction-induced cou- 
pling. Stroboscopic excitation "dresses" the atomic in- 
ternal states, so that each atom acquires a time averaged 
dipole moment of the order of d s = n 2 eaBf, oriented in 
the direction of E s . Even though the quantity / = At/T 
is assumed to be small, the induced dipole can be rather 
large for n ^S> 1. Taking for example At = Ins, T — 10/xs, 
and n = 20, we obtain d s = 0.1D. In the limit / — * 0, 
n 2 f = const, the resulting time dependent Hamiltonian 
can be replaced by its time average, leading to Eq.(jl]) 
with d = d s . A characteristic time scale in Eq.(jfl) is of 
order the inverse trap frequency u> . Hence, in our case 
the dynamics of the system is described by Eq.(|l]) with 
d = d S) if the condition At, T -C w^ 1 is satisfied. This 
has been tested numerically for At/T = 10~ 4 : We re- 
produced our previous results of the static GPE by solv- 
ing explicitly Eq.(|l|) for the stroboscopic "dressing" of 
atoms, i.e. by setting d — d s / f in the time intervals At 
and d = otherwise. 

Aside from inducing permanent electric dipoles, the 
stroboscopic dressing of atoms will somewhat modify the 
trapping potential and the scattering length related to 
the Van der Waals interatomic forces. The corresponding 
corrections will be proportional to the small parameter /. 
There will also be losses due to spontaneous emission and 
black body radiation |L4j . The rates of these processes 
become comparable with each other for n = 20 E3,E3, 
where the corresponding decay time is of order 20/is. In 
our scheme the lifetime will be thus ~ 20 /its// ~ 0.2 s. 

The laser resonant with a bare transition frequency 
"dresses" only the atoms that are sufficiently separated 
from their neighbors, since otherwise the dipole-dipole 
interaction shifts atomic resonances. Atomic pairs are 
"shielded" IT3] and not dressed at interatomic distances 
smaller than i?* , where the latter follows from the equa- 
tion HQ ~ d 2 R /R^, with fl = it/ At being the Rabi fre- 
quency associated with the "dressing" laser. For n = 20 
and At = Ins we have ft = 500MHz and R* = 0.7^m. 

The most dangerous "underwater stone" concerns in- 
elastic decay processes. Fortunately, the "shielding" can 
suppress Penning ionization: a strong suppression is ex- 
pected if atoms practically do not move during the short 
time At, and the distances at which the ionization oc- 
curs (~ rt 2 as) are significantly smaller than R*. For 
the parameters considered above this should be the case. 
The dipole-dipole interaction also induces the change of 



the effective Rabi frequency, and therefore the 2tt pulse 
condition is not strictly satisfied. Hence, a fraction of 
atoms remains in the Rydberg state between the stro- 
boscopic pulses and decays due to spontaneous emission. 
This fraction can be reduced by decreasing the quantity 
hd 2 At, where h is the gas density. A detailed analysis 
of inelastic processes in the conditions of stroboscopic 
dressing of atoms requires a separate investigation. 
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